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SUMMARY 


Tiis  report  presents  a  method  of  determining  a  confidence  integral  foe  a^«Ubi!ity  it  is 
esti gated  Uca  the  oean  tine  between  cqaipaecl  failures  and  the  mean  tine  to  repair  these 
failures. 

Then  the  failure  rate  aid  repair  rate  of  a  syslea  are  constants  independent  of  line,  it  can  be 
shown  that  the  estimates  of  mean  tine  between  failures  and  of  mean  tine  to  repair  the  failures 
can  be  combined  tc  have  tbe  variance  ratio  density.  The  estimated  availability  can  be  shown  to 
have  a  density  which  is  a  function  of  the  variance  ratio  density.  This  relation  can  then  be  used 
to  obtain  equations  foe  a  confidence  interval  for  availability. 

An  example  is  included  to  demonstrate  the  procedure  of  placing  a  confidence  interval  about 
the  estimated  avaiicbi'ity. 
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INTRODUCTION 


A;  the  U.  S.  Naval  !e  Center  (NtfC).  the  terns  “availability**  has  been  used  as  a  meas¬ 
ure  of  the  probability  that  a  weapon  sysiea  will  be  ready  at  a  prescribed  time  to  perform  its 
intended  Inaction.  In  this  sense,  availability  (A)  ls  defined  as 


0-6 


«*here  #  is  the  mean  time  between  e^eipsn^nt  failures  and  o  is  the  mean  time  to  repair  these  fail¬ 
ures.  The  actual  values  of  these  parameters,  though  unknown,  are  represented  hr  the  estimates 
A  A  A  t  A 

<?  and  o -  Since  6  and  &  are  calculated  from  data,  they  are  subject  to  random  variation  inherent 
in  sanding  The  estimate  of  availability 


A  A 
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is  tfces  also  subject  to  random  variation  and  can  be  expected  to  change  from  sample  tc  sample 
even  though  the  actual  value  of  availability  is  constant  over  these  samples. 

fn  the  past,  availability  has  been  reported  as  a  pom!  estimate,  os  single  value.  A  point 
estimate  of  a  parameter  is  net  very  meaningful  without  some  measure  of  the  possible  error  :n  the 
estimate.  An  estimate  of  a  paracseter  shorsld  be  accompanied  by  some  interval,  or  range  of  values, 
about  the  estir?&te  together  wath  some  measure  of  assurance  that  the  actual  value  of  the  parameter 
does  He  within  the  interval.  This  interval  is  commonly  known  as  a  confidence  interval. 

A  literature  search  o.d  not  disclose  any  documented  method  for  determining  a  confidence 

mtPival  for  availability.  T>e  purpose  of  this  repoit  is  to  document  such  a  method.  In  order  to  cal- 

rujate  a  confidence  interval  for  availability,  it  is  first  necessary  to  deteimine  the  probability 

A 

density  of  some  fraction  of  the  estimate  A.  This  report  is  specifically  concerned  with  deriving  3 

3  probability  density  function  under  some  limiting  assumptions  and  demonstrating  how  it  mav  fee 

A 

y*ed  to  calculate  the  limits  of  a  contidmce  interval  about  A. 

DERIVATION  OF  THE  DENSITY  FUNCTION  OF  0 

if  it  is  assumed  that  the  failure  rate  of  a  system.  -L  .  is  a  constant  independent  of  time,  then 
the  density  function  of  t.  (the  operating  time  between  the  (i-l)tn  and  the  ith  failure)  is 

f (t  )  -L  exp  I-  -±\ 

1  0  V\  0/ 

By  definition 


n 


where  n  is  the  rwr.ibe:  of  time  periods. 
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The  joint  density  function  of  (tj,  t2 . tR)  is 


By  transforoiation,  the  joint  density  function  of  (Q ,  t^t  t3,  .  .  tn)  is  obtained: 

V  *3 . *n>  =  a  (jY  C*P  (~  -f") 

A 

The  marginal  density  function,  g(Q).  is  obtained  by  integrating  over  the  variables  t^,  U . t^: 


(n-«5 


2 nn 

Let  the  random  variable  u  = -  .  By  a  change  of  variable,  the  density  function  of  g(u)  is 

obtained:  ^ 


2”  (n-1)  l 


'which  is  a  chi-squarcd  ( X 2)  density  function  with  2h  degrees  of  freedom. 


DERIVATION  OF  THE  DENSITY  FUNCTION  OF  6 

If  it  is  assumed  that  the  repair  rate  of  a  system,  -L-  ,  is  a  constant  independent  of  time,  then 

f5 

the  density  function  of  t.  (the  time  required  to  repair  the  jth  failure)  is 


By  definition 


j  =  l 


where  m  is  the  number  of  repairs. 

A 

The  density  function  of  <b  can  be  obtained  by  a  procedure  similar  to  the  one  used  for  the 
A 

dens!ty  function  of  0,  Thus 
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A 

g(6)  = 


(iJ 


tnm  exp 


(-4) 


(m-l) ! 


Let  the  random  variable  v  = 
obtained: 


g(v)  = 


1  e»P  (-  f) 


2a  (m-l) ! 


By  a  change  of  variable*  the  density  function  of  g(v)  is 


which  is  a  X2  density  function  with  2m  degrees  of  freedom. 


DEGREES  OP  FREEDOM 

In  the  derivation  of  the  density  function  of  the  estimate  0 ,  it  was  assumed  that  each  oper¬ 
ating  time  ti  terminated  in  a  failure.  Therefore*  there  were  n  values  of  t-  and  n  failures.  The 

A 

2n  Q  's 

random  variable  u  =  — 1—  is  distribute  as  XZ  with  2n  degrees  of  freedom. 

6 

In  the  event  that  the  last  operating  time,  tn*  does  not  terminate  in  a  failure  but  terminates  for 

some  other  reason,  such  as  the  end  of  a  monitoring  period,  and  0  is  calculated  as  the  total  oper- 

A 

ating  time  divided  by  number  of  failures,  (n-1),  then  the  quantity  u  =  2 (n-1)  6  ^  jjstfjbuted 

6 

as  X2  with  ?n  degrees  of  freedom. 

For  the  first  case,  the  number  of  degrees  of  freedom  associated  with  the  random  variable  u  is 
equal  to  twice  the  number  of  failures.  For  the  second  case,  the  number  of  degrees  of  freedom 
associated  with  the  random  variable  u  is  equal  to  twice  the  number  of  failures  plus  two. 

In  most  situations,  such  as  the  life  testing  of  equipment,  it  is  readily  apparent  which  one  of 
the  two  cases  applies.  If  the  test  is  terminated  at  a  predetermined  number  of  failures,  then  the 
first  case  applies.  If  the  test  is  terminated  at  a  predetermined  operating  time,  then  the  second 
case  applies. 

In  the  situation  where  monitoring  of  systems  :s  terminated  at  a  predetermined  calendar  time 
(for  instance,  one  week),  then  a  decision  must  be  made  as  to  which  of  the  two  cases  applies. 
Since  a  failure  occurs  at  a  discrete  point  in  time,  it  is  highly  unlikely  that  a  monitoring  period 
terminating  at  the  end  of  a  period  of  calendar  time  would  also  terminate  exactly  at  a  failure.  By 
Uus  argument,  the  second  case  would  apply. 


In  the  derivation  of  th^  density  function  of  .he  estimate  <£,  it  was  assumed  that  each  repair 

A 

time  t.  terminated  in  a  repair.  Therefore,  the  random  variable  v  -  _ _  is  distributed  as  X'4*  with 

J  d> 

2m  degrees  of  freedom. 
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CONFIDENCE  INTERVAL  ABOUT  A 


Since  u  and  v  are  independently  distributed  as  y2  with  2n  and  2m  degrees  of  freedom,  the 
quantity 

u/2n 

v/2m 

has  the  variance  ratio  density  (F)  with  2n  and  2m  degrees  of  freedom. 

The  following  relations  can  be  used  to  obtain  the  upper  (U)  and  lower  (L)  confidence  limits 

A 

of  a  confidence  interval  about  . 


CASE  1 


CASE  2 
A 

/0_\  =  _n_  j£_  p 

W  /y  n-1  A  ;  2n,  2m 


A 

/  <f>  \  _  n  d>  p 

\~6  )L  n-1  A  -J-  ;  2n,  2m 


where 

n  =  the  number  of  failures  in  case  1 
n  -=  the  number  of  failures  plus  one  in  case  2 
m  -  the  number  of  repairs 
l  -  a  -  the  level  of  confidence 

f  -  the  upper  -ll5L  fractile  of  the  cumulative  F  density  with  2n  and  2m 

A—  ;  2n,  2m  2 

2 

degrees  of  freedom 

F  -  the  lower  —  fractile  of  the  cumulative  F  density  with  2n  and  2m 

V  ;  2n.  2m  2 

degrees  of  freedom 
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If  availability  is  expressed  as 


then  the  upper  and  lower  limits  of  a  confidence  interval  about  A  may  be  obtained  from  the 
following  relations: 


EXAMPLE  OF  DETERMINATION  OF  A  CONFIDENCE  INTERVAL 

A  typical  situation  at  NMC  may  involve  the  monitoring  of  an  F-4B  squadron  during  a  one- 
week  weapons  training  exercise.  The  following  is  an  example  of  how  a  confidence  interval  may 
be  determined  from  data  collected  during  one  of  these  exercises.  Given  the  following  fictitious 
data: 


Total  equipment  “on”  time . 200  hours 

Total  number  of  failures . 14 

Total  time  expended  in  repair . 100  hours 

Total  number  of  repairs . 20 


then 


0  =  .200  =  14.286 
14 


$  =  i00  =  5.000 

20 


A  =  — —  =  0.741 
A  A 

0  4-0 

For  a  95  per  cent  confidence  interval  on  availability 


(11)  (AM.)  (1.94)  -  0.727S 

\  14  /  \  14.286/ 


*It  is  possible  for  there  to  be  more  repairs  than  failures  if  several  aircraft  begin  the  exercise  with  failure 


A 

n  ^  p 

n-1  A  0  025  <30’  40> 

0 


1 14  A  14.286 

/  \  2.01 

1 

1 

1 

^ _ J 

4 

1.1806 

1 

1 

-  0.8428 


l  J±\  17275 


=  0.5789 


